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ABSTRACT 

Optical microcavities support Whispering Gallery Modes (WGMs) with a very high quality factor Q. 
However, WGMs typically display a far-field isotropic emission profile and modifying this far-field profile 
without spoiling the associated high Q remains a challenge. Using a 2D annular cavity, we present a 
procedure capable to achieve these two apparently conflicting goals. With the correspondence between 
the classical and the wave picture, properties of the classical phase space shed some light on the charac- 
teristics of the wave dynamics. Specifically, the annular cavity has a well separated mixed phase space, 
a characteristic that proves to be of crucial importance in the emission properties of WGMs. While the 
onset of directionality in the far-field may be achieved through parametric deformation pQ of the distance 
cavity-hole centers, d , this contribution presents a method to control the emission profile via a second 
parameter, the hole radius rrj. The influence of the classical dynamics to control and predict the field 

T 1 emission will be demonstrated. 
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O 1. INTRODUCTION 

Whispering Gallery Modes (WGMs) are known to exist in certain types of microcavities. For instance, 
the class of regularly shaped dielectric structures includes the sphere, the toroid and the disk (thin 
cylinder) cavities. The strong confinement of the light field, characteristic of the WGMs, appears as 
sharp peaks of the spectral response of these dielectric microcavities. Accordingly, these cavities possess 
high quality (high Q) modes. Although this property is of central importance for most, if not all, 
applications of microcavities, the complete uniformity of the emitted far-field distribution of WGMs 
I remains an impediment for useful light source applications (e.g. microlasers [21 [3]). 

^. This paper is the sequel of a previous contribution (see [I] in these Proceedings) . We will be interested 

in the parametrical control of the far-field patterns of a WGM located inside an annular dielectric cavity. 
Our approach follows the observation made in [Tj that the behavior of the far-field of a WGM may be 
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triggered from uniformity to non-uniformity using the center-to-center distance d (see figure |l(a) I . We 



o 

will exploit the correspondence between the rise of non-uniformity in the far-field distribution (wave 
picture) and the extent of non-regular dynamics in classical phase space (ray picture) to predict and 

o 



control the emission. 

2. CLASSICAL CONCEPTS AND GENERAL CONSIDERATIONS 



The geometry of the annular dielectric cavity [4j is depicted in figure 1(a) It consists of an outer disk 



of radius R a and refractive index n^ in which a circular inclusion of radius r and refractive index n^ is 
embedded and whose center is a distance d from the center of the main disk. The cavity is surrounded 
by a dielectric medium of index n . In the classical limit, the dynamics within the cavity amounts to 
follow individual rays as they propagate freely in between collisions with the outer boundary and the 
inclusion. For an ensemble of many trajectories, the phase space is constructed from the successive 



coordinates {si,pi} (figure 1(b)) recorded on the surface of the outermost boundary: {si} are the arc 
length and {pi — sin^i} where Xi is the angle incidence with respect to the normal. Note that the 
rays are specularly reflected at the outermost boundary whereas we allow transmission (|p| < nh/nc) 



or reflection (\p\ > nh/nc) at the inclusion interfaces. The result is displayed in figure (1(c)). Several 
features are worth mentionning, some of them unique to the annular cavity. 
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The system under study possess a mixed dynamics. Even more, phase space is well divided into two 
principal regions: a regular domain, \p\ > pnr = (d+ ro)/Ro , filled with trajectories, reminiscent of the 
integrable disk dynamics, oblivious to the inclusion, and a non-regular domain, \p\ < pnr where irregular, 
mostly chaotic, trajectories are observed. The extent of the region where the dynamics is non-regular is 
simply governed by the position d and size rg of the inclusion and proportional to d+rg. Another domain of 
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Figure 1: Geometry and phase space of the annular cavity, (a) A circular inclusion of refractive 
index n^ and radius r$ is placed a distance d from the geometrical center of the cavity. The distance of 
the inclusion Rq — (d + ro) to the boundary will be kept constant in the following; (b) Classical phase 
space coordinates, arc length s (alternatively angular position cj> — s/ Rq) and p — sin^; (c) Phase space 
constructed from the impact coordinates {4>i,Pi\ for different initial conditions. Two distinct regions, con- 
taining regular and irregular trajectories respectively, are apparent and well separated by the geometrical 
condition \p\ = Pt$r = (d + tq)/ Rq. The boundary of a TIR region (\p\ > ptir ~ 0.3) is also indicated. 



interest is the classical emission (or escape) region defined by the E = {(s,p) : < s < 2irR , \p\ < Ptir}, 
under the Snell-Descartes condition of Total Internal Refexion (TIR) at the outer boundary 



Ptir = n /n c 



(1) 



Clearly for j>nr > Ptir, the classical far-field emission is strictly determined by the non-regular "chaotic" 
trajectories pQ. Moreover, only a small subset of E is actually accessible for emission, we call this subset 
the effective emission region, and denoted it by W. It consists of the trajectories initially located outside 
of E that get mapped to E through a single iteration of the (implicit) Poincare map, V : (si,Pi) i— > 
(si+i, Pi+i) ■ The region W is then defined as 

W = V{V~ x {E)nE) (2) 

with the complement E = {(s,p) : < s < 2ttRq, ptir < \p\ < !}■ Trajectories entering E are thereafter 
bound to escape the cavity. The dynamical system so-defined does not conserve energy (i.e. the total 
amount of initial conditions (or initial intensity) decreases as the trajectories reach E). 

The domain W is obviously modified by parametric variations of both d and ro. Since we have observed 
that the far-field emission of a WGM becomes non-uniform with increasing d (for fixed ro) [3[S], and 
accordingly increasing pnr, we propose that the directionality of the emission may be further engineered 
by keeping d + r^ (pnr) constant and changing ro, thereby modifying the "mixing" properties of the non- 
regular region without affecting its size. The quantitative modifications induced will now be investigated 
in the wave dynamics under the guidance of the previous classical considerations. 

3. EFFECT OF PARAMETER r ON THE FAR-FIELD BEHAVIOR 

We approximate our physical annular cavity by a planar (2D) system and search for resonant modes 
(here TM polarization) by solving Hclmholtz equation 



[V 2 +n 2 (r»fc 2 ]^(r,^)=0 



(3) 



under appropriate boundary conditions. The solutions are obtained by a careful implementation of the 
scattering formalism of [B]. The physical parameters of the system are set as nc = 3.2, rih = n = 1, 
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Figure 2: Far-field intensities, full curves: combined envelope of the two symmetries of mode (11, 1), 
IV'fn i)| 2 + iV'ni nl / dashed curves: non-coherent sum of the classical trajectories reaching the far-field 
for 3 values of ro, (a) r$ = 0.064i?o, (b) ro = 0.127i?o„ (c) r$ = 0.211i?o, while d + r$ is kept fixed 
at 0.55i?o. For ease of comparison, the intensities in each panel have been arbitrarily normalized to the 
maximum of the classical intensities. 



Ro = 1 keeping the control group d + ro constant at = 0.55 Ro. The important regions are then delimited 
by Ptir = 0.3125 and pnr = 0.55. For the sake of demonstration, we have selected the high Q WGM 
(11,1) at kR sa 4.5. The notation (to, n) is meant to identify the angular momentum number and 
the number of radial nodes respectively, and k is the wavenumber. The particular choice of d + ro is 
motivated by the high far-field contrast value C\\ for ro = 0.20i?o and d — 0.35i?o OP- This contrast 
measure C m „ varies between and 1: meaning that the observed far- field has fully retained its original 
uniform mo character, and 1 indicating that the perturbed mode has mixed with other angular momentum 
components and that its field distribution is highly non-uniform. We have retained the two symmetries 
of mode (11, 1) in the full-wave calculations (even/odd, VVii i-p relative to the center-to-center axis). 

We follow the mode (11, 1) as ro is varied from 0.05i?o to 0.40i?o and although its quality factor drops 
from about 10 to roughly 10 5 over the interval, it remains the dominant resonance, albeit with a changing 
character. The far-field emission is recorded as a function of the observation angle 9 from the non-coherent 
sum of \ip? n x J 2 and \tp? n jJ 2 - The far-field intensities are displayed in figure (|2l for 3 separate values 
of ro. The emission profiles evolve from a sharp contribution near 9 = it at ro = 0.064i?o (figure 2(a) ) 
to an increasingly broad multi-peaks envelope with equally important contributions around 9 = 7r/4, w 
and 37r/4 at ro = 0.127i?o (figure 2(b)) and to a depleted emission in the backward direction in favor 
of two important forward peaks at ro = 0.211i?o (figure 2(c)). In short, the objective of modifying the 
directionality of the far-field emission of a WGM through parametric control is achieved. 

To connect with the intuition gained in the previous section, we have carried out a ray-escape simulation 
[7] in the classical representation. Assuming that the transition probability to the effective escape region 
W is highest near the j»nr limit, initial conditions are chosen uniformly in a thin strip of phase space with 
V ~ Pnr- For a given ro, each initial condition (so,po) is followed through consecutive impacts (si,Pi) 



with the outer boundary where a reflected intensity Zj+i = Zj(l 
transmission coefficient, 

T(P) = 



\T{P; 



i+U 



v/i^r 



y/l - p 2 + ^/{n /n c ) 2 -p 2 



is calculated with Fresnel 



(4) 



Iterations stop when the intensity associated with a trajectory has dropped from its initial value Io to 
an arbitrarily small number. The escaping intensities of every trajectory are then binned according to 
observation angles 9 to form the far-field distribution. Finally, the classical far-field is convoluted with a 
gaussian function of width a = 27r/4mo to smooth out irregularities. The width is chosen to account for 
the angular resolution of the pure WGM (too, n) which has 2to lobes in the azimuthal direction. 

The comparison with the full-wave far-field profiles is presented in figure [2] At this low value of 
nckRo ~ 15, it is no surprise that the agreement is only qualitative although the main features are 
present in both calculations. Our exploratory simulations are nevertheless encouraging and we expect 



the procedure to deliver better correspondence at higher values of nckRo. One should stress however 
that the classical approach is universal and generic. In the annular geometry studied here, a large number 
of high Q WGM (m, n) are available above \p\ = pnr and they are well localized around their resonance 
condition pwgm = rn/ (nckRo). As long as j>wgm > Pnr, they are candidate for far-field directionality 
while preserving their near-field uniformity, i.e. low-loss, high Q character, and phase space design is a 
possible route to optimization. 

Figure ([3]) gives a sense of the range of values of the control parameters where one can modify (engineer) 
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Figure 3: Constrat measure C\\ versus ro. For fixed d+ro = 0.55 (distances in unit ofRo), 2 regions, 
non-uniform and uniform far-field, arbitrarily divided at C\\ = 0.5, define the range of ro for which the 
far-field profile has a certain character. The 3 dots correspond to the same values of ro presented in figure 
(2). 

at will the far-field distribution keeping a high quality factor. The contrast measure Cn(ro) stays above 
0.80, meaning strong mixing with other angular components m ^ m , for r 6 [~ 0, ~ 0.25]. As r is 
further increased, Cn drops further until it reaches low values where the WGM has recovered its regular 
uniform character. 



4. CONCLUSION 

In this paper we have shown that it is possible to control the directionality of the far-field emission 
through the modification of the radius ro of a circular inclusion in a annular dielectric cavity. Evidences 
have been presented that the far-field profile of both full-wave and classical simulations show similar 
structures. Furthermore, for a given annular geometry, the classical far-field profile is universal. Our 
results along with those of [11 [5] pave the way to design scenarios for high Q directional emission from 
WGMs of disk-shaped microcavities. 
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